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Abstract: Masu salmon, Oncorhynchus masou masou, is one of the most valuable fishery species that has been introduced to China, 
though to date no studies on the genetic diversity and genetic relationship among hatchery populations has been performed with 
molecular markers. We undertook such a study and sampled 120 individuals from three hatchery stocks and analyzed 20 
microsatellite loci. All loci were polymorphic and a total of 91 alleles were detected. A relatively low level of genetic diversity was 
revealed with effective number of allele of 3.1094, 3.3299 and 3.1894 and expected heterozygosity of 0.6600, 0.6648 and 0.6638 in 
the three stocks, respectively. Deviations from Hardy-Weinberg equilibrium were found due to heterozygote deficit. Accordingly, 
evidence of genetic bottlenecks were found in the three stocks. An individual assignment test demonstrated that 85% of individuals 
were correctly assigned into their original stocks. Pairwise Fst revealed that significant differentiation occurred between these three 
stocks. The results of the study indicated that disequilibrium of genetic structure and differentiation has occurred in all three stocks. 
This information collectively provides a basis for measures to avoid of loss of genetic diversity and introgression in Chinese 


aquaculture. 
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Masu salmon, O. masou masou, were widely 
distributed in Japanese mountain streams (Kano et al, 
2010) before being introduced to China in the 1990s 
(Wang, 1998). The fish is now one of the most valuable 
fishery species in cold-water aquaculture. At present, 
commercial catches of O. masou masou in Japan have 
decreased and hatchery programs have been established 
to release fry to its habitat (Yu et al, 2010). In China, 
because of its good taste and flesh quality, the fish has 
been an important aquaculture species (Wang, 1998). 
However, successive artificial selection and serious 
stress from inbreeding have led to the remarkable 
degradation in the fish: early sex maturity, smaller size 
and silvered body color (Zhang et al, 2011). A second 
batch of O. masou masou was introduced into China in 
2008 to counter these changes. 

High genetic diversity is important for the population 
conservation because it is closely related to population's 
fitness and adaptability to environment as well as 
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selection (Filgueira et al, 2010). Loss of genetic variation 
among hatchery stocks is a common phenomenon caused 
by genetic drift in many species, including mrigal (Aung 
et al, 2010), Japanese flounder (Sekino et al, 2002) and 
common carp (Kohlmann et al, 2005). At present, there 
are three O. masou masou hatchery stocks in China, but 
their genetic diversity and relationship have not been 
reported despite clear indications of genetic 
abnormalities that prompted the second importation of 
stock in 2008. To avoid a repeat of the earlier incidence, 
information on genetic variation of these hatchery stocks 
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is urgently needed to develop a coherent management 
and conservation plan for this species. 

Microsatellites DNA markers have been proven to 
be useful markers in researching genetic diversity and 
population structure in fish such as Japanese flounder 
(Sekino et al, 2002) and common carp (Kohlmann et al, 
2005). We aim to investigate the genetic diversity and 
population relationships of the Chinese three stocks of O. 
masou masou using 20 pairs of microsatellite primers 
developed for rainbow trout to assess the genetic 
diversity and their relationship within and among stocks 
of O. masou masou. 


MATERIALS AND METHODS 


Sample collection and history of stocks studied 

Two batches of fertilized eggs of O. masou masou 
were introduced into China, the first in 1996 and the 
second 2008. Both were kept in Bohai Experimental 
Station. In 1998, parts of O. masou masou individuals 
were taken to Beijing Shuntong Rainbow Trout Culture 
Center and regarded as a separate stock population after 
several generations of artificial propagation, resulting in 
three stocks. We named the three stocks we used 
BO1996, Beijing, and BO2008. In 2010, we collected 40 
samples from each of the presumptive stocks. And the 
samples of the BO2008 stock were collected directly 
from the introduced population. Fin-clip tissues were 
preserved in 75% ethanol (Jia et al, 2008) for BO1996 
and Beijing stocks. BO2008 stock samples were 
collected as fry just before beginning to float. 


DNA extraction and microsatellite DNA analysis 

Genome DNA was extracted using the standard 
proteinase K, phenol-chloroform-isoamyl procedure 
(Wei et al, 2001). 20 microsatellite loci (Tablel) 
developed for other salmon were used in the study. The 
OMM3005 primer was from NCBI 
(http://www.ncbi.nlm.nih.gov/nuccore/20142355) ^ and 
others were from the articles of Rexroad et al 
(2002a,b;2005). Primer conditions, mainly annealing 
temperatures, were optimized for amplification. PCR 
was performed using a reaction volume of 15 uL 
containing 50 ng of genomic DNA, 0.2 umol/L of each 
dNTP (TaKaRa, China), 0.1 umol/L of each primer, 1.5 
mmol/L MgCh, 1x PCR buffer (TaKaRa, China) and 1U 
Taq DNA polymerase. The thermal cycling (Biometra) 
began with 3 min at 94 ?C, followed by 27 cycles of 
denaturing at 94 ?C for 30 s, annealing temperature from 
50 °C to 62 °C depending on different primers for 30 s, 
elongation for 30 s at 72 °C, followed by a final 
extension for 5 min at 72 ?C. PCR products were 
separated in 10% polyacrylamide gels. 
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Statistical analyses 

The effective number of alleles per locus (Ne), 
number of detected alleles (Na), observed heterozygosity 
(Ho), expected heterozygosity (He) and allele frequency 
were calculated using POPGENE 1.32(Yeh et al, 2000). 
To test the difference between populations of Ne and He, 
non-parametric t-tests (Sokal & Rohlf, 1995) were 
performed using Microsoft Excel. 

Deviations from Hardy-Weinberg | equilibrium 
(HWE) were estimated using GENEPOP 4.0 (Raymond 
& Rousset, 1995). Exact P-values were evaluated with 
Markov chain method (dememorization=10000; 
batches-100; iterations per batch-5 000). Additionally, 
deviations from HWE in the direction of heterozygote 
excess or deficit were also carried out. 

Analysis of molecular variance (AMOVA), 
implemented in GenAIEx (Peakall & Smouse, 2006), 
was used to partition total genetic variation hierarchically 
among stocks. 

Genetic distance between populations was estimated 
by the PHYLIP software package, version 3.69 
(Felsenstein, 1993). Bootstrap re-sampling (1 000 
replicates) was performed to test the reliability of 
dendrogram implemented in PHYLIP software. 

STRUCTURE 2.2 (Pritchard et al, 2000) was used 
to determine the most likely number of genetic clusters 
(K=1-5). An assignment test was also performed in 
order to determine the extent to which individuals could 
be correctly assigned to their population of origin. In the 
analysis, the simulation was used assuming admixture 
and correlated allele frequencies between populations 
with a 10 000 replications burn-in period and 100 000 
MCMC replicates. 

Evidence of population bottlenecks were tested with 
Bottleneck 1.2.02 (Cornuet & Luikart, 1997) and the 
Wilcoxon sign-rank test was used to determine if a 
population exhibits a significant number of loci with 
heterozygote excess with a two-phased model of 
mutation (TPM, as recommended in the program).The 
mode-shift indicator was also employed to discriminate 
bottlenecked populations from stable populations. 


RESULTS 


PCR results of 20 pairs of microsattelite primers 
According to the clarity and size of banding 
patterns, 20 pairs of primers showed good images and all 
were polymorphic in the sampled populations (Figure 1). 


Genetic diversity and linkage disequilibrium 
Summary statistics showing allelic variation within 
loci and population are presented in Table 1. A total of 91 
alleles were detected in the 20 loci among 120 
individuals with an average of 4.55, ranging from three 


Wwww.zoores.ac.cn 


Genetic diversity and differentiation of masu salmon between and within cultured populations inferred from microsatellite DNA analysis — E35 


b bb bb bbbbb* *, bt wibi nb, bbbp^bp?* p 


ceddiede,.aqd^dcCdd d dd,ddedeed de^ 





eece 


BO1996 


27 UL. 


b 
€ ddddd 


aa bbb 
dddddd°° ddd dddddddd 


um. b 





BO2008 


b 
did cdd 


' dd d,ddd d d 
c, aa 


> f 





Figure 1 Polyacrylamide gel results at the OMM5017 locus for three stocks of masu salmon 


a, b, c, d and e means different alleles. 


alleles at loci AF352740, AF352758 and OMMS5020 to 7 
alleles at locus AF352763. All loci were polymorphic in 
all samples from three stocks. Ho and He within stocks 
ranged from 0.2973 (BO2008) to 0.9714 (BO1996) 
and from 0.3225 (Beijing) to 0.8290 (Beijing), 
respectively. 

Exact tests for genotypic linkage disequilibrium 
showed significant deviations in the stocks of BO1996 
(OMM1036 and OMM1035, OMM1036 and AF352744, 
AF352738 and AF352750 pairs) and Beijing (OMM1045 
and AF352763, AF352763 and af352746, AF352750 and 
AF352744 pairs) at the 5% level. But disequilibria were 
not consistently observed in the same locus pairs in each 
population. 


Genetic structure with and among populations 

The effective number of allele per population was 
3.1094, 3.3299 and 3.1894 in BO1996, BO2008 and 
Beijing respectively. The mean expected heterozygosity 
per population was 0.6172, 0.6171 and 0.6299 in 
BO1996, BO2008 and Beijing respectively. No 
significant difference of genetic diversity was found 
between stocks. 

Deviations from HWE were found in all three 
populations (Table 1). Heterozygote deficit was observed 
in BO2008 stock while for the other two, however, 
heterozygote deficits existed at only few loci and did not 
achieve population level. The BO1996 population 
exhibited the highest number of loci (11) deviated from 
HWE and then BO2008 (10) and Beijing (9) stock 
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populations. Two loci (OMMS5017 and OMMS010) 
deviated from HWE in all populations but a heterozygote 
deficiency existed in different population (OMMS017 for 
BO2008 and Beijing, OMMS010 for BO2008 at 5% 
level). A significant bottleneck was found for three 
stocks (P [heterozygote excess] — 0.00000, 0.00001 and 
0.00136 for BO1996, BO2008 and Beijing, respectively 
and the shifted model). 

For the real number of clusters, once the cluster (K) 
is reached, L(K) at larger Ks plateaus or continues 
increasing slightly. In our analysis, L(K) at larger Ks 
plateaus when K=3, which is consistent with the groups 
as individual sampled (Figure 2). And the 4K value that 
could give a clear peak at the true value of K according 
to procedures of Evanno et al (2005) was in Figure 3. 
The self-assignment test resulted in correct assignment of 
102 of the 120 individuals in the three sampled groups, 
corresponding to 85.0% (Figure 4). The number of 
misclassified individuals in the three groups was four, 
seven, and seven with the ratios of 3.33%, 5.83%, and, 
5.83% respectively. 

The locality of origin for each individual is 
indicated in the X-axis. 


Genetic relationships 

The global test for population differentiation was 
significant with Fst value of 0.055 (P=0.01). Pairwise 
Fst estimates and level of significance for all population 
pairs are presented in Table 3. All Fst values were 
significant (all P=0.01) and the highest Fst was found 
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Table 1 Genetic variability of three O. masou masou stocks at 20 microsatellite loci 


Locus Repeat Population P Ne Ho He 

BO1996 0.0013* 3.9801 0.7500 0.7582 

OMMU L (GT)17 BO2008 0.0077* 3.8186 0.6250 0.7475 
BEIJING 0.0317* 2.4078 0.6000 0.5921 

BO1996 0.0000* 3.9851 0.9000 0.7585 

SNIN (GT)25 BO2008 0.0000* 4.3185 0.5250 0.7782 
BEIJING 0.0017* 3.8415 0.8250 0.7491 

BO1996 0.0988 3.3058 0.8250 0.7063 

OMM5003 (GA)44 BO2008 0.0146* 2.4298 0.7500 0.5959 
BEIJING 0.4824 2.7997 0.7250 0.6509 

BO1996 0.0084* 2.9331 0.6750 0.6674 

OMM3008 (GT)19 BO2008 0.2133 2.4600 0.4865 0.6016 
BEIJING 0.0058* 2.8855 0.8500 0.6617 

BO1996 0.2681 3.0149 0.8205 0.6770 

OMM5900 (CA)18 BO2008 0.0003* 2.5499 0.8462 0.6157 
BEIJING 0.0475* 3.3934 0.8250 0.7142 

BO1996 0.0174* 3.7559 0.4750 0.7430 

MMI? (TCTA)12 BO2008 0.0003* 2.3145 0.2973 0.5757 
BEIJING 0.2520 3.0274 0.5500 0.6782 

BO1996 0.0111* 3.7736 0.7000 0.7443 

QMM1044 (TCCA)14 BO2008 0.0000* 5.2605 0.5526 0.8207 
BEIJING 0.1396 4.5403 0.6154 0.7899 

BO1996 0.9844 3.1232 0.7692 0.6886 

OMMI036 (TATC)18 BO2008 0.5301 4.8784 0.7368 0.8056 
BEIJING 0.1108 3.5242 0.8750 0.7253 

BO1996 0.0620 4.5134 0.6923 0.7885 

OMMI (TAGA)31 BO2008 0.0508 5.1200 0.6500 0.8149 
BEIJING 0.0000* 4.3478 0.9500 0.7797 

BO1996 0.5363 3.5455 0.6923 0.7273 

DIST (CA)14 BO2008 0.0000* 3.9661 0.4872 0.7576 
BEIJING 0.0388* 3.3465 0.5789 0.7105 

BO1996 0.3406 1.8827 0.6053 0.4751 

Me (CT)6N2(CT)4 BO2008 0.1241 2.2467 0.7179 0.5621 
BEIJING 0.3765 2.2651 0.6667 0.5658 

BO1996 0.0424* 3.5794 0.6250 0.7297 

OMMa003 (CA)37 BO2008 0.3607 4.3770 0.6667 0.7816 
BEIJING 0.0479* 4.5007 0.8500 0.7877 

BO1996 0.7945 1.7847 0.4250 0.4453 

AF3239 (TAGA)15 BO2008 0.3727 2.5457 0.7500 0.6149 
BEIJING 0.3994 1.9606 0.5526 0.4965 

BO1996 0.1491 4.3986 0.7714 0.7839 

AF352763 (GATA)26 BO2008 0.0579 6.3908 0.7436 0.8545 
BEIJING 0.0802 5.4688 0.9714 0.8290 

BO1996 0.1019 2.3529 0.5250 0.5823 

AF332130 (GATA)14 BO2008 0.0926 2.4317 0.6923 0.5964 
BEIJING 0.0726 2.4634 0.7250 0.6016 

BO1996 0.0115* 2.8594 0.5000 0.6589 

AF332146 (TATC)7 BO2008 0.0070* 2.9304 0.4750 0.6671 
BEIJING 0.1266 3.3508 0.6000 0.7104 

BO1996 0.0112* 2.8889 0.4615 0.6623 

RENS (GATA)19 BO2008 0.0000* 2.7112 0.3077 0.6394 
BEIJING 0.0558 3.0239 0.5128 0.6780 

BO1996 0.0000* 2.6644 0.8500 0.6326 

eo (AC)14 BO2008 0.2775 2.2028 0.5897 0.5531 
BEIJING 0.0000* 2.8803 0.7750 0.6611 

BO1996 0.0114* 1.9987 0.7179 0.5062 

AT332 158 (GATA)13 BO2008 1.0000 1.8376 0.4865 0.4621 
BEIJING 0.4266 1.4672 0.3000 0.3225 

BO1996 0.9234 1.8470 0.5128 0.4645 

DEUS (CA)17 BO2008 0.0457* 1.8075 0.6154 0.4525 
BEIJING 0.0000* 2.2923 0.9500 0.5709 

BOI996 | — -—-— 3.1094 0.6647 0.6600 

Mean .  --—— BO2008 |  -—— 3.3299 0.6001 0.6648 
BEIJING  -—- 3.1894 0.7149 0.6638 


P: estimate of probability of deviation from Hardy-Weinberg equilibrium; *: significant departure from Hardy-Weinberg equilibrium (P<0.05); Ne: the effective 
number of alleles; Ho: observed heterozygosity; He: expected heterozygosity. 
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between the BO1996 and BO2008 populations. Global 
subdivision with the use of AMOVA revealed that the 
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percentage of variation among populations was 5%, 4% 
among individuals, and 91% within individuals. This 
indicates that the total variation of microsatellite loci was 
due to variation among individuals, rather than among 
stocks. 
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Figure 4 Proportional membership (Q) of each individual in the three clusters 


Table3 Pairwise Population Fst Values at 20 
microsatellite loci 


Pop ID BO1996 BO2008 Beijing 
BO1996 gie 0.01 0.01 
BO2008 0.062 ad 0.01 
Beijing 0.050 0.052 gie 


Fst Values below diagonal. Probability values based on 99 permutations are 


shown above diagonal. 


Pairwise genetic distances among hatchery 
populations (Table 4) ranged from 0.1282 (BO1996— 
Beijing) to 0.1349 (BO1996-BO2008), which was 
consistent with Fst results. 


Table 4 Nei's genetic distance (below diagonal) 


Pop ID BO1996 BO2008 Beijing 

BO1996 LM RR ZZ 

BO2008 0.1671 n 

Beijing 0.1282 0.1349 seek 
DISCUSSIONS 


Our study provides useful insights into the genetic 
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diversity of O. masou masou and the data indicates 
relative low levels of genetic diversity but significant 
population differentiation the three O. masou masou 
stocks. 


Within-population diversity and population structure 
Low levels of allele diversity among O. masou 
masou were suggested by the 20 microsatellite loci. The 
number of alleles ranged from 3.1094 to 3.3299 per locus, 
compared to the average value from 12 species of 
freshwater fish (4—9.1) (De Woody & Avise, 2000) and 
from O. masou masou populations in Japanese and 
Korean populations (A=6.42—9.71) (Yu et al, 2010). 
However, the observed and expected heterozygosity 
(Ho=0.6001-0.7149 , He-0.6600—0.6648) was 
moderate, slightly higher than the average of freshwater 
fish (H=0.54 averaged of 13 species) (De Woody & 
Avise, 2000) and lower than O. masou masou 
(Ho=0.698-0.817, He-0.722-0.814) (Yu et al, 2010). 
For three stocks, samples were introduced as fry and they 
may have comprised a few full-sib families rather than 
representatives of their original populations (Na-Nakorn 
et al, 2004), which resulted in the low level of alleles. 
We also used primers developed for other salmon. Cross- 
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amplified primers and the loss of low-frequency alleles 
after several generations of artificial propagation may 
result in low genetic variation. BO1996 and Beijing 
population showed even lower genetic diversity 
compared to BO2008 in terms of number of alleles and 
heterozygosity. BO1996 and Beijing populations 
experienced a longer period of being cultured since 
introduction into China as compared to the BO2008 
population. This time difference might lead to a loss of 
diversity as a result of inbreeding and genetic drift. 

In general, a small number of alleles are a signature 
of population bottleneck (Norris et al, 1999) due to allele 
loss. The present results indicate that significant genetic 
bottlenecks were detected with the shifted model in all 
populations, which implies that a dramatic loss of alleles 
happened in our samples. In the next breeding or 
hatcheries, a sufficient. number of individuals are 
necessary to avoid serious bottleneck. 

Deviations from HWE were found in three 
populations due to heterozygote deficiencies, likely caused 
by null alleles, a mixture of independent populations, 
non-random mating, or artificial or natural selection. In 
this study, half of loci deviated from HWE in the 
sampled populations and heterozygote deficits did not 
occur in all the loci and populations, which did not 
reflect features of null allele effect but likely of mixture 
of independent populations (for BO2008 stock), non- 
random mating, or artificial or natural selection (for 
BO1996 and Beijing stock). 


Genetic differentiation between stock populations 

Genetic differentiation was observed between/ 
among different stocks. In all inter-stock tests, deep 
differentiation was found between BO1996 and BO2008 
and then Beijing and BO2008, which was consistent with 
the different stocks’ origin. 

Significant population differentiation was observed 
between BO2008 and the other two. As a newly 
introduced stock, BO2008 might be genetically different 
from original populations. It is interesting that two other 
stocks, BO1996 and Beijing, are also differentiated 
significantly from each other, as Beijing samples were 
recruited from BO1996 stocks. This could be due to 
genetic drift phenomenon such as non-random mating, 
inbreeding or selection that occurred in hatcheries. 
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